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We present and investigate an analogue model for a controllable photon generation via the dy¬ 
namical Casimir effect (DCE) in a cavity containing a degenerate optical amplifier (OPA) which 
is pumped by an amplitude-modulated field. The time modulation of the pump field in the model 
OPA system is equivalent to a periodic modulation of the cavity length, which is responsible for 
the generation of the Casimir radiation. By taking into account the rapidly oscillating terms of 
the modulation frequency, the effects of the corresponding counter-rotating terms (CRTs) on the 
analogue Casimir radiation emerge clearly. We find that the mean number of generated photons and 
their quantum statistical properties exhibit oscillatory behaviours, which are controllable through 
the modulation frequency as an external control parameter. We also recognize a new phenomenon, 
the so-called “Anti-DCE ”, in which pair photons can be coherently annihilated due to the time- 
modulated pumping. We show that the Casimir radiation exhibits quadrature squeezing, photon 
bunching and super-Poissonian statistics which are controllable by modulation frequency. We also 
calculate the power spectrum of the intracavity light field. We find that the appearance of the side 
bands in the spectrum is due to the presence of the CRTs. 
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I. INTRODUCTION 

One of the most important and interesting predictions 
of quantum field theory is that the quantum vacuum state 
is not truly empty, but instead filled with virtual parti¬ 
cles which are continuously created and annihilated due 
to the zero-point fluctuations [1]. In quantum electrody¬ 
namics these fluctuations lead to a number of measurable 
effects in the realm of microscopic physics, such as the 
van der Waals force between atoms in vacuum, natural 
widths of spectral lines. Lamb shift, and anomalous mag¬ 
netic moment of electron [1]. In macroscopic level, the 
vacuum fluctuations manifest also themselves in a variety 
of observable mechanical phenomena and the archetype 
of them is the so-called static Casimir effect [2, 3]. This 
effect which arises from the pressure that virtual photons 
exert on stationary boundaries, has been observed in dif¬ 
ferent experiments [4, 5] and plays a very important role 
both in fundamental physics investigations as well as in 
understanding the basic limits of nanomechanical tech¬ 
nologies [6, 7]. 

Another manifestation of the zero-point fluctuations is 
the creation of real particles out of the vacuum induced 
by the interaction with dynamical external constraints 
[8]. In this effect, energy is transferred from the ex¬ 
ternal field to the zero-point fluctuations transforming 
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them into real particles. Although a material boundary 
is somewhat analogous to an external constraint, there is 
no particle creation from the vacuum in the case of static 
boundaries (static external constraints). However, if the 
boundaries are nonstationary, and the boundary condi¬ 
tions depend on time, there is particle creation from the 
vacuum in addition to the Casimir force. Such quan¬ 
tum vacuum amplification effect is commonly referred 
to as the dynamical Casimir effect (DCE) [9, 10] (for an 
extensive recent review, see [11]). Whereas the static 
Casimir effect arises from a mismatch of vacuum modes 
in space, the DCE originates from a mismatch of vacuum 
modes in time domain. Examples of the DCE range from 
cosmology, such as particle creation in a time-dependent 
cosmological background [12-14] to cavity QED, such as 
photon generation in Fabry-Perot cavities with moving 
mirrors [15-21]. There is an extensive literature on var¬ 
ious aspects of the creation of photons due to the DCE 
using various theoretical methods. For example, pho¬ 
ton creation was described based on a Hamiltonian ap¬ 
proach [22]. Cavities with the insertion of dispersive 
mirrors or a slab with a time-dependent dielectric per¬ 
mittivity were considered [23-26]. Multiple-scale anal¬ 
ysis was applied to the calculation of the flux of cre¬ 
ated particles [18]. Photon creation in a harmonically 
oscillating one-dimensional cavity with mixed boundary 
conditions has been analyzed [27]. The DCE has also 
been investigated in the context of cavity QED [28, 29], 
in Bose-Einstein condensates [30], in excition-polariton 
condensates [31], in superconducting circuits[32], and in 
a quantum-well assisted optomechanical cavity [33]. Be- 
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sides its intrinsic importance as a direct proof of the vac¬ 
uum fluctuations, various practical applications of the 
DCE has been proposed, e.g., for generation of pho¬ 
tons with nonclassical properties [34-37], for high preci¬ 
sion optical interferometry [38], for generation of atomic 
squeezed states [39], for multipartite entanglement gen¬ 
eration in cavity networks [40], and for quantum commu¬ 
nication protocols [41]. 

From the experimental point of view, the detection of 
the DCE induced by moving boundaries is a serious prob¬ 
lem. The main difficulty is largely due to the challeng¬ 
ing prerequisite of generating appreciable mechanical fre¬ 
quencies of moving boundaries (of the order of GHz) to 
obtain a detectable number of photons. Although there 
are some experimental proposals for the detection of the 
DCE involving real mechanical motion of boundaries (see 
for instance [42]), several analogue systems have been 
proposed for observing the Casimir radiation based on 
simulating moving boundaries by considering material 
bodies with time-dependent electromagnetic properties. 
The first analogue system consisting of a nonlinear op¬ 
tical medium with time-dependent refractive index was 
introduced in [9] and later on, its experimental imple¬ 
mentation was reported [43]. In [44] the Casimir radi¬ 
ation within a cavity containing a thin semiconducting 
film with time-dependent conductivity and centered at 
the middle of the cavity has been investigated. An ex¬ 
plicit analogy between photon emission of a coherently 
pumped thin nonlinear crystal of type inside a cavity 
and the photon generation via the DCE have been estab¬ 
lished [45]. In [46] an experimental layout for measuring 
the DCE has been proposed where the optical length of a 
cavity is varied in time by means of a train of ultrashort 
pulses that modulate the effective refractive index of the 
cavity material. By using a varying boundary condition 
in a superconducting waveguide [47, 48] which simulates 
a single moving mirror, the first experimental demonstra¬ 
tion of the DCE has recently been reported [49]. More¬ 
over, a second observation of the DCE has been reported 
[50] based on periodical changes in the index of refraction 
of a microwave cavity with a Josephson metamaterials. 

Motivated by the above-mentioned investigations on 
analogue models for the DCE, in the present contribu¬ 
tion we propose a model system for the generation of 
the Casimir radiation in a cavity containing a degenerate 
optical parametric amplifier (OPA) which is pumped by 
an amplitude-modulated field according to the harmonic 
law 


ep{t) = e + r]cosnt, (1) 

where £ is a constant amplitude and rj, H are, respec¬ 
tively, the amplitude and frequency of modulation. The 
temporal change in pump amplitude leads to a modula¬ 
tion of the cavity length which is equivalent to an ap¬ 
parent displacement of the cavity mirrors. In all of the 
previous investigations carried out on the DCE, a kind 
of rotating wave approximation (RWA) has been applied 
to neglect the contribution of rapidly oscillating terms 


OPA 



K 


Spit) 





7 

>■ 

(b) 

J 




Mechanical 

analogue 


i5L(t) = 2Lo^si7ifH 



FIG. 1. (Color online) (a) Schematic representation of a 
Fabry-Perot cavity containing an OPA which is pumped by 
an amplitude-modulated field, (b) The analogy between the 
cavity with time modulated pumping and the cavity with os¬ 
cillating mirror 


of the modulation frequency( see for instance [51-55]). 
The main feature of our treatment is that we release this 
limitation to include the corresponding counter rotating 
terms (CRTs). We will find that taking into account the 
CRTs in the model under consideration leads to an in¬ 
terference between the constant amplitude, £, and the 
oscillatory term. The total number of generated photons 
consists of three contributions, i.e., constant pumping, 
time varying pumping, and their interference which orig¬ 
inates from CRTs. Due to this interference, the mean 
number of generated photons and their quantum statis¬ 
tical properties exhibit oscillatory behaviours which are 
controllable through the modulation frequency Q. 

The paper is structured as follows. In section. II we 
describe the model system and show the analogy between 
the time modulation of the pump and the time modula¬ 
tion of the cavity length in the conventional DCE. We 
also solve analytically the quantum Langevin equations 
for the cavity-field operators. We calculate the mean 
number of generated Casimir photons and analyse its 
temporal behaviour in section. III. We study the quadra¬ 
ture squeezing, photon counting statistics and the power 
spectrum of the generated Casimir radiation in section. 
IV. Finally, we summarize the main results of the paper 
in section. V. 


II. DESCRIPTION OF THE SYSTEM 

As schematically shown in Fig. 1(a), we consider a sys¬ 
tem composed of a sub-threshold degenerate OPA inside 
a Fabry-Perot cavity with fixed mirrors at zero temper¬ 
ature. We assume that the OPA is pumed by a coupling 
field whose amplitude is time-modulated according to Eq. 
(1). We also take into account the dissipation of the cav¬ 
ity field to the environment with a decay rate k. We 
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show that this non-stationary cavity QED results in the 
DCE. In a degenerate OPA a pump photon of frequency 
u}p is down- converted into a pair of identical photons of 
frequency wq. The Hamiltonian describing the system in 
the interaction picture is given by 


Hi = ^pit) (a 1'2 - a^) , ( 2 ) 


where a{a^) is the annihilation (creation) operator of the 
cavity field. Here, it should be emphasized that in the 
above Hamiltonian we will not make use of the RWA 
for the term cosHt in the amplitude ep{ t) and we want 
to examine the influence of the corresponding CRTs on 
the generation and quantum statistical properties of the 
Casimir photons. Before proceeding to discuss the dy¬ 
namics of the system under consideration, it is worth¬ 
while to beriefly describe the analogy between the time 
modulation of the pump and the time modulation of the 
cavity length in the conventional DCE. Following the 
same method as used in [45], it is easy to show that the 
time-modulated nonlinear parametric interaction in the 
model system, described by the Hamiltonian (2), is equiv¬ 
alent to an apparent periodic displacement of the cavity 
mirrors which is given by SL{t) = 2LQ{r]/n) sin fit, where 
the effective length of cavity is Lq = Lcav + K'^OPa ~ 1 )j 
with I and no pa as the length and the refractive index 
of the nonlinear crystal, respectively (see Fig. 1(b)). Fur¬ 
thermore, the maximum of the effective speed of the mir¬ 


ror is obtained as (u/c)max 


5L{t)/c 


T)—Kj2 


kLq jc (for 


example, if Lq ^ cm and n ~ lO"* —10® Hz, then {vlc)max 
is of order 10“"* — 10“^). As we will show latter, in com¬ 
parison to the analogue model introduced in [45], where 
the refractive index of nonlinear crystal coated on the 
mirror is modulated, our analogue model has the advan¬ 
tage of providing the possibility of controlling the genera¬ 
tion of the Casimir photons and their quantum statistical 
properties by adjusting the frequency of modulation, D. 
This controllability results from taking into account the 
CRTs in the amplitude ep{t). 

The quantum Langevin equations of the system, taking 
into account the cavity dissipation, can be written as 


^ = -^a + ep{t)d^ +Fc{t), 
^ = -fat+e,(t)« + Ft(i), 


where Fc{t) is the noise operator associated with the 
dissipation of the cavity field with zero mean value, 
i.e.,(^Fc{t)'j = 0. For a cavity mode damped by a vacuum 

reservoir, the noise operator in zero temperature satisfies 
the Markovian correlation functions: 

= (4) 


(Ft(t)F,(t')) = (F,(t)F,(t')) = {FKt)FKt')) = 0 . 

„ . „ ^ 

By defining the operators X := + d^Y := d^ — d, and 

F± := F^{t)ztFc{t) the quantum Langevin equations can 


be written as 

1= (ep(t)-f)A + F+(t), 

Y = -{epit) + ^)Y + FUt). 

The solution of the set of Eqs. ( 6 ) is given by 

Xit) = m (X(0) + f* dt'F+{t')g,{t')) , 

Yit) = / 2 (f) (y(0) + f* dt'F.{t')g 2 {t')) , 

where 

= g(e-K/ 2 )t-|-r)sin(nt)^ gi(f) = 

/ 2 (f) = g-(e+«/ 2 )t-'^sin(at)^ 52 (f) = . 

( 8 ) 

with 77 = ry/D. Thus, the held operators are given by the 
following expressions 

a(f) = i(/i(f)l(0)-/2(f)y(0) + 

/i(f) /g <if'F+(f')5i(f') - / 2 (f) /o dt'F_{t')g 2 {t')^ , 
d\t) = \[h{t)X{Q) + h{t)Y{Q)+ 

hit) /J dt"F+{t'')ght") + hit) Jo df"F-(f")52(f")) • 

(9) 


III. MEAN NUMBER OF GENERATED 
PHOTONS 


We are now in a position to calculate the mean num¬ 
ber of intracavity photons by using the analytic solution 
given by Eq.(9). Assuming that the cavity held is ini¬ 
tially prepared in the vacuum state,| 0 ), the mean number 
of photons at time f is given by 


(h(f)) = 0.25 


/iW (l + ^fodt'ght')^ 
+/ 2 W (l + fifodt'glit')^ 



( 10 ) 


In order to get a clear insight into Eq. (10), we rewrite 
it as follows 

(nit)) = 0.25 [(e-'>'+*-^“"(^*) -P e-J-t+zsin(nt) _ 2 ) 

+K (/+(f) -I- /-(f))], 

( 11 ) 

where z = 2fj, 7 ± = k ± 2 e, and /±(f) = 

/±(f)e-'^±*T2zsin(nt) 


I±(t)=Io 


gUn 

gfit') 


dt' = 


E 


7 j- -\-imVL 


o{l±-\-'irnVL)t _ 


1) 


( 12 ) 


in which Jm is the Bessel function of the hrst kind of 
integral order m. By using the properties of the Bessel 
function[56] we obtain 
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/±(i) = u;± Jo(* 2 :) (e t'±* - 1 ) + Y. [2mn sm{2mnt) + j± {cos{2mV,t) - e '>'±‘)] 

^ ^ (13) 

T E -yl+(Ztli)'in 2 [7± Sin((2m + l)nt) + (2m + l)n (e '^±* - cos((2m + 1)^))] 


where w± = 1 /"f±. In this manner, the mean number of 
photons can be separated into three contributions as 


(7r(t)) — ^OPA “b Prj “b Pinterf erence — '^OPA PC< 

where 

^ --Y+i ^ t 

noPA = -+ —e - -—e , 

7+7- 27 + 27_ 


astmir 5 


(14) 

(15) 


Tin = - 1 ) 


1 1 

-1- 

7+ 7- 


p-7+t p-l-t 


7+ 


7- 


(16) 


4n. 


interference = - l) (1 - g-T'-*)) + {e^- l) (e-^+‘ + (1 - 6-^'+*)) 


-l-fc V ( p+2sin(nt) ) 2J2rri(i^) 

oo ^ , , 


pmH sin( 2 mllt) + 7 ± {cos{2mVtt,) — e '•'='=*)] 


(17) 


T E nl+{ 2 m-i)^a 2 b± sin((2m - l)Vlt) + (2m - 1)0 (e - cos((2m - 1)^))] 


As can be seen, riinterference is an oscillatory function 
of time and it can be shown that in the sum over m 
in Eq. (17) the term corresponding to m = 1 has the 
dominant contribution. We remind that since OPA is 
operating below threshold, the parameters e and rj are 
constrained by the inequality k > 2(e + 7 ). In the ab¬ 
sence of the time modulation, i.e., 7 = 0 , both the terms 
Tin and riinterf erence Vanish and the mean number of gen¬ 
erated photon reduces to nopA given by Eq.(15). In 
the steady state, itopa reaches (2 £^)/(k^ — 4£^). Eig- 
ure 2 (a) illustrates no pa versus the dimensionless time 
Kt for different values of e/n. As is seen, with increas¬ 
ing e/k towards its maximum value(< 0.5) the rate of 
photon generation increases. In Eigs. 2(b)-(d) we have 
plotted ncasimir = Pr, + Pinterf erence aS a funCtion of the 
dimensionless time fit for different values of tt/fl (with 
a fixed value of k). As can be seen, ncastmir exhibits 
an oscillatory behaviour in the course of time. The os¬ 
cillations can be attributed to the CRTs which manifest 
themselves as an interference between the constant and 
time varying parts of the external pumping Ep. Eor a 
given value of e/ 7 , decreasing the modulation frequency 
ll(or equivalently, increasing n/Q. for a fixed value of n) 
increases the amplitude of oscillations. In addition, it is 
seen that for a fixed value of k/O the mean number of 
generated Casimir photons increases with increasing e/ 7 . 
Figures 2(b)-(d) also show that, depending on the values 
of /c/11 and e/ 7 , the mean number of Casimir photons 
takes negative values within some time intervals, which 
means that (h(t)) < no pa, although we have always 
{n{t)) > 0. This effect can be described as a coherent 
annihilation of pair photons due to the time modulation 


of the external pumping. With increasing the value of O, 
the mean number of annihilated photons decreases and 
for a fixed value of fl this annihilation occurs for large 
values of e/ 7 . The effect of the photon creation and anni¬ 
hilation can be attributed to the exchange of the energy 
between the external modulated pump and the quantum 
vacuum fluctuations. If the energy is transferred to the 
quantum vacuum fluctuations, photon pairs are created. 
Conversely, if the energy is gained from the quantnm 
vacuum, photon pairs are annihilated. Indeed, the term 
cos(llt) corresponds to the energy exchange such that the 
external pump varies between e — 7 and e -I- 7 . It is inter¬ 
esting to note that this annihilation effect has also been 
reported very recently [25] in a microscopic toy model 
for cavity DCE consisting of a non-stationary dielectric 
slab inside a fixed single-mode cavity. The slab, which is 
modelled as a set of N non-interacting atoms, oscillates 
according to an external law, while its dielectric prop¬ 
erties are modulated externally via electric or magnetic 
fields. It has been shown that for = 1, the modula¬ 
tion of atomic parameters can lead to coherent annihi¬ 
lation of three photons accompanied by the creation of 
one atomic excitation. The authors has called this effect 
“Anti-DCE”. 


An advantage of the system discussed so far, in com¬ 
parison to the proposed system in [45], is the controlla¬ 
bility of the mean number of generated Casimir photons 
through the external parameters fl, e, and 7 . 
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FIG. 2. (Color online) (a) Mean number of generated pho¬ 
tons in the absence of time modulation {rj = 0 ) versus the 
dimensionless time nt for different values of e/k. (b), (c) and 
(d) Mean number of Casimir photons versus the dimension¬ 
less time at, for different values of n/fl and e/rj. Here, we 
have set Ac/27r = 10‘^Hz. 


IV. QUANTUM STATISTICAL PROPERTIES 
AND POWER SPECTRUM OF GENERATED 
PHOTONS 

In this section, we study the quantum statistical prop¬ 
erties of the generated photons, including quadrature 
squeezing and counting statistics. We also calculate the 
power spectrum. In particular, we will show how one 
can manipulate these properties through controlling the 
parameters of the external modulated pumping. 


A. Quadrature Squeezing 

The squeezing property of the radiation field can be 
analysed by calculating the variances of the quadrature 
operators. Using Eq.(9) and properties of the noise op¬ 
erator Fc{t), given by Eqs.(4) and (5), the variances of 
the quadrature operators di and 02 are, respectively, ob¬ 
tained as 

(Adi(t))^ = (1 + kJ* dt'gf(t')) , 

(Aa2(t) f = -b k/p dt'gl{t')j , 

where di = X/2 and 0,2 = lY/2 are quadrature opera¬ 
tors and obey the commutation relation [oi, 02 ] = 1 / 2 . A 
state of the radiation field is said to be squeezed when¬ 
ever {Adi{t))'^ < 0.25, {i = 1, 2). The degree of squeezing 
can be measured by the squeezing parameter Si (z = 1 , 
2), defined by Si{t) = A{Adi{t))^ — 1. Then, the condi¬ 
tion for squeezing in the quadrature component can be 


simply written as Si{t) < 0. Using Eq.(12) the different 
contributions in Eq.(18) can be separated in a similar 
manner as the case of photon number in the previous 
section. In Eig. 3(a) we have plotted the squeezing pa¬ 
rameters Si and S 2 versus the dimensionless time Kt for 
the case of the constant pumping (77 = 0 ) and for two 
values of e/n. As is seen, the quadrature component 02 
exhibits squeezing at all times, while the component di 
is never squeezed. Furthermore, as time goes on, the 
squeezing parameters Si and S 2 tend to asymptomatic 
values 2 £/(k — 2e) and —2e/{K + 2e), respectively. Nu¬ 
merical results corresponding to the temporal evolution 
of the squeezing parameters Si and S 2 in the case of 
time-modulated pumping (77 ^ 0 ), for different values of 
k/Q, e/k and 77 /K (k being constant) are shown in Figs. 
3(b) and 3(c), respectively. The effect of CRTs leads to 
the oscillatory behaviour of quadrature squeezing in the 
course of time such that the dominant frequency of these 
oscillations depends on Q. Unlike the constant pump- 
ing(Fig. 3(a)), in this case both quadratures di and 02 
can be squeezed in some time intervals (see Figs. 3(b), 
(c)). It should be reminded that this effect is due to 
the CRTs. As shown in Figs. 3(b) and(c), the quadrature 
squeezing can be controlled by the modulation frequency, 
Q, as an external control parameter. For fixed values of k, 
rj, and e when Q decreases {k/Q increases) the maximum 
degree of squeezing and the amplitude of oscillations in¬ 
crease. Furthermore, for a fixed value of k/Q, when e/rj 
increases, the degree of squeezing of 02 increases while 
the squeezing of the component di is suppressed. Our 
calculations show that if k/Q be very smaller than unity, 
the temporal behaviours of the quadrature variances are 
similar to those obtained for the OPA but with very small 
oscillations. By controlling the modulation frequency Q 
the maximum quadrature squeezing is about 0.13, i.e., 
50% noise reduction is achievable. We remind that we 
have plotted the squeezing parameters for the total gen¬ 
erated photons. However, it is easy to show that the 
quadrature variances for generated Casimir photons are 


Adi ^Casimir) — (Aui) (Au^^g) 


(Au^,/;) T (Ah2,mier/erence) ; 


(19) 


where {Ad^^e)'^, {Adi^^Y, and {AdiwinterferenceY are, re¬ 
spectively, the contribution of constant pumping, the 
time modulated term and the interference term which 
is due to the CRTs. Numerical results reveal (not shown 
here) that the behaviour of {Adi^casimirY' is almost simi¬ 
lar to that of {AdiY, fiiat is the difference is quantitative 
rather than qualitative. 


B. PHOTON COUNTING STATISTICS: 
Autocorrelation function and the Mandel parameter 

The autocorrelation function is proportional to the 
probability of detecting one photon at time t t given 
that another photon was detected at earlier time t. It is 
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FIG. 3. (Color online) (a) The time evolution of the squeezing parameters Si{t) and 52 (t) versus the dimensionless time Kt in 
the case of r; = 0 (constant pumping) and for two values of e/fc. (b) The time evolution of the squeezing parameter Si{t) versus 
the dimensionless time Qt for the case of ry 7 ^ 0 (time-modulated pumping) and for different values of k/H, e/k, and rj/n. (c) 
The time evolution of the squeezing parameter 52 (t) versus the dimensionless time Qt for the case of t; 7 ^ 0 (time-modulated 
pumping) and for different values of k/Q, e/k, and rj/K. Here, we have set k/2'k = 10^ Hz. 


defined by 

( 2 ) ^ {a^t)d\t + T)a{t + T)a{t)) 

When the radiation field satisfies the inequality (t) < 
g^^^(O), the photons tend to distribute themselves prefer¬ 
entially in bunches rather than at random(photon bunch¬ 
ing) . On the other hand, if g^^^ (r) > g^^^ (0) fewer photon 
pairs are detected close together than further apart (pho¬ 
ton antibunching). Using the Gaussian properties of the 
noise forces and the quantum regression theorem [58, 59] 


we get, after algebraic manipulations, the expression of 
autocorrelation function 


= 1 + 


(a'l'(t)d'l' (t -I- t)) 

(d^t)a(t)f 


+ 


(^a(t)^d(t + t)^ 
(dHt)d(t)f 


( 21 ) 

In order to determine two correlation functions appear¬ 
ing in the right hand side of Eq. (21) one can use Eq. 
(9) together with the correlation functions for the noise 
operators Fc, given by Eq. (4) and (5). The results are 
as follows 


(dl (t)d(t + r)) 


0.25 fi(t)fi(t + T) (^l + Kfgdt'g^(t')J + f2(t)f2(t + t) (i + k/o dt'ffi(t')) 
_2e-0-5'‘l^lCosh (s 


1 £T -I- 2 ^Sin^Cos 


( 22 ) 


(d^t)d^t + T)) =0.25 fi(t)fi(t + T) + + f 2 (t)f 2 (t + r) (^k dt'gl(t') 

_2g-o.5«;|r| 2 ^ sin(l^) cos( ^^^2^^^ 


- 1 


(23) 


Considering these two expressions for times much longer 
than the cavity relaxation time (t ^ k~^) we have plot¬ 
ted ill Fig- 'l(a) for the case of constant 

pumping (g = 0 ) as a function of kt, for two values of 
s/k (k = constant). As is seen, in this case, the gener¬ 
ated photons exhibit bunching phenomenon. Also, de¬ 
creasing e leads to an inhibition of photon bunching. In 
Fig. 4(b), we have plotted foi' l^e case 

of time-modulated pumping (g 7 ^ 0 ) as a function of 
Qt for different values k/Q and e/g. As is seen, for a 
fixed value of k/Q when e/g increases photon bunching 
decreases and g^‘^\T)/g^'^\Q) approaches 0.3. Further¬ 
more, with increasing k/Q (decreasing the modulation 


frequency 12 ) when k is constant, photon bunching in¬ 
creases and g^'^^r )/approaches 0. The fast decay 
of the autocorrelation function is due to the CRTs. Nu¬ 
merical analysis shows that If k/Q < 1, the behaviour 
of the autocorrelation function is similar to that ob¬ 
tained for the OPA, but it has oscillatory behaviour at 
long delay times. Increasing s/g leads to the decreas¬ 
ing of the amplitude of these oscillations. The manifes¬ 
tation of photon bunching is attributed to the possibil¬ 
ity of emission of two photons simultaneously by DCE 
that is an OPA like process (it can be described by 
the Hamiltonian of Eq.(2)). Up to now, we have found 
that for a fixed value of k with the decreasing value 
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of r2, the mean number of generated Casimir photons, 
the maximum degree of squeezing and photon bunch¬ 
ing increase. To determine the photon statistics, one 
can calculate the Mandel parameter which is defined as 
Q = ((An)^ — (h))/ (h) = ((a’^^a^) — (h)^)/ (h) where 

8 

^ (24) 

Z = 1 


the functions Fi{t) and Ai{t) are as follows 
-I- 3 (/i^(t) + f2'^{t)) 

A2{t) = A^{t) = AQ{t) = A'j{t) = i — /2^(t)), 

A4{t) = l{h^{t) + h\t)-2e-^), 

As{t) = 1, 

(25) 


Fi{t) = l, 
F2{t)=Fr{t) = 
F^{t)=Fi{t) = 


f (/i (0 /o 9l{x)dx - fi{t) J* gl{x)dx^ , 
F^it) = FQ{t) = I nfUt) /o glix) 

+«/! W Io9i2ix)dx - 2(1 - e"”*) 


-Fs)^) = TV.P2 + ^ 6^3 + F5F4. 


(26) 

For Q > 0 (Q < 0); the statistics is super-Poissonian 
(sub-Poissonian); Q = 0 stands for Poissonian statis¬ 
tics. Figure 5(a) shows the temporal evolution of the 
Mandel parameter in the absence of time modulation 
(77 = 0 ) versus the scaled time Kt for different values 
of e/k. As is seen, as time goes on, the Mandel param¬ 
eter increases and it is finally stabilized at an asymp¬ 
totic value. With the increasing value of e the super- 
Poissonian statistics is enhanced. Furthermore, the rates 
with which the reaching the asymptotic values occur are 
inversely proportional to s; the smaller e is, the more 
rapidly the Mandel parameter tends to the asymptotic 
value. Figures 5(b) and (c) show the temporal evolu¬ 
tion of the Mandel parameter for the generated Casimir 
photons versus the scaled time fit for different values of 
k/Q, e/k and rj/n. As is seen, for k/H > 1, Q{t) > 0 
for all times, which means that the generated Casimir 
radiation obeys super-Poissonian statistics. The Mandel 
parameter exhibits an oscillatory behaviour due to the 
CRTs such that the frequency of oscillations depends on 
n. Furthermore, we observe that the super-Poissonian 
behaviour can be controlled by SI and e/tj in such a 
way that if k and SI are constant, by increasing e/ry one 
not only decreases the amplitude of the oscillations of 
the Mandel parameter, but also decreases its maximum 
value. Furthermore for a fixed value of e/ry < 1, in¬ 
creasing of /c/S1(k =constant) leads to decreasing of the 
amplitude of the oscillations in the Mandel parameter 
and its maximum value(Fig.5(c)). On the other hand, 
for e/yy > 1, increasing of k/S1(k =constant) leads to in¬ 
creasing of the amplitude of the oscillations in the Man- 
del parameter and its maximum value(Fig.5(b)). For 
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FIG. 4. (Color online) (a) Normalized autocorrelation func¬ 
tion , versus the scaled delay time kt in the 

absence of time modulation and for different values of e/k. (b) 
Normalized autocorrelation function , g^^\T)/(0), versus 
the scaled delay time Qt in the presence of time-modulated 
pumping and for different values of k/Q, and e/g. Here, we 
have set k/2tt = 10"^ Hz. 


K/n ^ 1 , the Mandel parameter has not a clear oscil¬ 
latory behaviour. In this case, as is seen, at the very 
initial stages of evolution {nt <C 1), Q{t) < 0, which 
means that the generated Casimir radiation obeys sub- 
Poissonian statistics but then it increases towards posi¬ 
tive values (super-Poissonian). 


C. Power Spectrum 

We next proceed to calculate the power spectrum of 
the intracavity light held. The power spectrum can be 
expressed in terms of the Fourier transform of the two- 
time correlation function in steady state as follows 

Sioj) = -Re dTF-^{dHt)ait + T))^^, (27) 

Jo 

where the two-time correlation function,(d^(t)a(t-I- t)) 

, is given by Eq.(22). Figures 6(a)-(f) illustrate the 
normalized power spectra for different values of e, Q, rj 
and 9s = fits (scaled interaction time in the steady 
state) and constant value of K/27r = lO^Hz. In Fig. 
6(a) we have plotted the normalized power spectrum 
S' 7 v(w) = S'(a;)/S'(e) versus H/k for ry = 0 and for two 
values of e. As is seen, increasing e not only causes the 
strength of the peak of the spectrum to increase, but 
also leads to narrowing of the spectrum. In hg.6(b)-(f) 
we have plotted the spectrum of the generated photons 
inside the cavity in the presence of time modulation as 
a function of w/fl. As can be seen, for hxed value of 
k/H, the width of the spectrum decreases and its height 
increases if e/ry increases. In general, for constant value 
of K the width of the spectrum increases if k/VL increases 
or the modulation frequency decreases. With varying 
the scaled time, 9s, the shape of the spectrum, i.e., its 
width, height and number of peaks, vary as a result of 
CRTs that cause the spectrum to be a function of time ts- 
Furthermore, varying 9s leads to the appearance of the 
symmetric side peaks in the spectrum with same width 
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FIG. 5. (Color online) (a) Time evolution of the Mandel parameter, Q{t), as a function of the dimensionless time Kt, in the 
absence of time modulation for different values of e/ft. (b) Time evolution of the Mandel parameter, Q{t), as a function of the 
dimensionless time Qt, in the presence of time modulation for e/n = 0.3 and tj/k = 0.1 and different values of k/Q. (c) Time 
evolution of the Mandel parameter, Q{t), as a function of the dimensionless time Qt, in the presence of time modulation for 
e/fv = 0.1 and tj/k = 0.3 and for different values of k/Q. Here, we have set k/2tt — lO'^ Hz. 
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FIG. 6. (Color online) (a) Normalized power spectrum, 5']v(tu) = S(tj)/S'(e), in the absence of time modulation, i.e., ?7 = 0, as 
a function of oj/k for different values of e/k. (b)-(f) Normalized power spectrum, S'jv(u;) = S{u})/S{Q) versus uj/Q for different 
values of e, t], Q, and 6s = Qts. Here, we have set k/2tt — 10“^ Hz. 


and height which depend on the cavity decay rate k, con¬ 
stant amplitude e, modulation amplitude r] and modula¬ 
tion frequency Ql. The distance between side peaks and 
their locations depend on the modulation frequency. The 
appearance of the side bands is due to the interference 
terms corresponding to the CRTs. Therefore, the spec¬ 
trum can be controlled by adjusting the modulation fre¬ 
quency and Os- 

It may be interesting to compare the intracavity spec¬ 
trum obtained above with the spectrum of the output 
photonic field. In fact the input-output theory[60] links 


the emitted output field a-out to the intracavity field a by 
a simple relation 


aout{t) = - a^n{t), (28) 


where, a{t) is the intracavity light field operator and 
din(t) = (l/v^)^c(t) is the input field. After algebraic 
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manipulations, we find field operators without using the RWA for rapidly oscil¬ 

lating terms of the modulation frequency, we have found 
+ t) \ = K {a^{t)a{t + t)) - {t)F^{t + t)\ that the mean number of intracavity photons consists 

' l~. \ \ 'of three contributions: constant pumping, time varying 

pumping, and their interference resulting from the CRTs. 


(29) 

At zero temperature, all of the correlation functions that 
appear in Eq. (29) are zero except for the first term, and 
thus 

{a\^^(t)aout(t + = K {a^{t)a{t + t)) . (30) 

This implies that the spectrum of the output field is 
proportional to the spectrum of the intracavity light 
field, i.e., Souti^) = kS{uj). Therefore, the Casimir ra¬ 
diation can be detected from the comparison between 
the output field spectrum in the absence of the time 
modulation(77 = 0) and the spectrum in the presence of 
the time modulation. 


V. SUMMARY AND CONCLUSIONS 

In summary, we have proposed an analogue model 
for controlled generation of photons via the DCE in a 
nonlinear cavity containing a degenerate OPA which is 
pumped by a time-modulated field. In the system un¬ 
der consideration the time modulated nonlinear para¬ 
metric interaction is analogous to an apparent periodic 
displacement of the cavity mirrors. By solving analyti¬ 
cally the quantum Langevin equations for the intracavity 


The results reveal that the temporal behaviour of the 
mean number of generated photons can be controlled via 
the parameters of the time-modulated pump field. We 
have also identified the so-called “Anti-DCE”, in which 
pair photons are coherently annihilated due to the time- 
modulated pumping. 


In addition to the generation of Casimir photons, 
we have discussed analytically and numerically their 
quantum statistical properties including the quadrature 
squeezing and photon-counting statistics. We have also 
examined the power spectrum of the generated Casimir 
radiation. We have found that the generated radiation 
exhibits quadrature squeezing, photon bunching, and su¬ 
per -Poissonian statistics in the course of time, which 
are controllable through the parameters of the time- 
modulated pumping. Furthermore, the main charac¬ 
teristics of the power spectrum, i.e., its height, width, 
and number of peaks can be controlled, in particular, by 
adjusting the modulation frequency and the interaction 
time in the steady state. The results also reveal that the 
comparison of the output power spectrum in the presence 
of the time- modulated pumping with one which obtained 
by constant pumping provides a way for the detection of 
the generated Casimir photons. 
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